Scrambling and Gate-Induced Fluctuations in Realistic Quantum Dots 
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We evaluate the magnitude of two important mesoscopic effects using a realistic model of typical quantum 
dots. "Scrambling" and "gate effect" are defined as the change in the single-particle spectrum due to added 
electrons or gate-induced shape deformation, respectively. These two effects are investigated systematically in 
both the self-consistent Kohn-Sham (KS) theory and a Fermi liquid-like Strutinsky approach. We find that the 
genuine scrambling effect is small because the potential here is smooth. In the KS theory, a key point is the 
implicit inclusion of residual interactions in the spectrum; these dominate and make scrambling appear larger. 
Finally, the gate effect is comparable in the two cases and, while small, is able to cause gate-induced spin 
transitions. 
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I. INTRODUCTION 

An important way to characterize quantum dots (QDs)!^ 2 ^ 
the simplest artificial nano-structure with electrons quantized 
in all three dimensions, is by the parametric evolution of their 
properties. The most common external parameter is mag- 
netic field because of its flexibility of tuning, 1 but other pa- 
rameters are also used. Here we are concerned with the ef- 
fect of changing the electron number 7Y or the external gate 
voltage U g , referred to as the scrambling and gate effects, re- 
spectively, in Coulomb blockade (CB) experiments 
The most striking feature of the CB regime is sharp peaks 
in the conductance through the quantum dot as a function of 
gate voltage. As shown in Fig.^ at each conductance peak, 
the number of electrons residing in the dot changes by one; 
across a peak spacing, the gate voltage changes to bring an- 
other electron into the dot, and deforms the confining poten- 
tial in the meantime. The scrambling and gate-induced shape 
deformation effects were both introducedii*^ in connection 
with experiments on the spacing between CB conductance 
peaks, 4 - 6 ' 7 ' 8 ' 9 and have also been used to interpret CB peak 
height correlationspiS 

The scrambling and gate effects can both be quantified 
through the variation in the single-particle spectrum of the 
system, {£;}. Since electron-electron interactions are impor- 
tant for quantum dots in the Coulomb blockade regime, one 
must clearly consider the effect of such interactions on the 
single-particle spectrum. Here we evaluate the scrambling and 
gate effects using both density functional theory and Thomas- 
Fermi calculations for realistic geometries of quantum dots. 
We address two main issues: 

First, while the magnitude of the scrambling and gate ef- 
fects has been estimated for hard-wall quantum dots coupled 
to a large gate (one which deforms the entire dot), 1314 exper- 
imental quantum dots have, of course, smooth confining po- 
tentials, and are typically deformed with a narrow "plunger" 
gate. We evaluate these experimental features using our real- 
istic model of quantum dots, showing that they influence the 
magnitude of the scrambling and gate effects strongly. 

Second, what "single-particle spectrum" should one use in 
evaluating the scrambling and gate effects? Roughly there 



are two types of single-particle spectra that can be defined in 
an interacting system. The first is a spectrum from a self- 
consistent mean field theory such as Hartree-Fock (HF)A£ or 
Kohn-Sham spin-density functional theory (KS-SDFT)ii&The 
second is the spectrum of a reference Hamiltonian which con- 
tains the interactions only at a smooth (classical-like) level. 
The most natural choice is the eigen-spectrum of the effec- 
tive potential calculated from Thomas-Fermi (TF) theory; this 
constitutes the Strutinsky approach (S-TF)piii^The difference 
between these two types of spectra is familiar, for instance, 
from discussions of the meaning of the eigenvalues in the 
self-consistent approach: Recall that the self-consistent eigen- 
value is related to the energy for removing an electron from 
that level and that a sum over such eigenvalues double counts 
the interaction energy among those electrons, 1619 neither of 
which is true for the eigen-spectrum in the reference Hamilto- 
nian approach. 

The meaning and magnitude of the scrambling and gate 
effects depend on which type of single-particle spectrum is 
used. We emphasize that this is not a question of which ap- 
proach is the more accurate, but rather of what part of the 
fluctuations of the total energy is assigned to these effects. 
For instance, when using a reference potential as in S-TF, the 
fluctuations as a parameter changes associated with interac- 
tions are separated into two distinct parts. The first comes 
from changes in the single-particle energies as the smooth 
Thomas-Fermi potential varies - we consider these the "gen- 
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FIG. 1 : (Color online) Illustration of scrambling and gate effects in 
Coulomb blockade conductance fluctuations. 
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uine" scrambling and gate effects. The second contribution 
involves the screened Coulomb interaction between meso- 
scopic density fluctuations. This "residual interaction" term 
corresponds to the weak interaction between Landau quasi- 
particles in a Fermi liquid picture. In a self-consistent ap- 
proach such as KS-SDFT, they appear in the self-consistent 
energies ef S (N); thus parametric evolution of the KS levels 
involves both genuine scrambling and residual interactions. 

The paper is organized as followings. In the next section, 
we outline briefly the two theoretical approaches, i.e. the 
Kohn-Sham method and Fermi-liquid like Strutinsky approx- 
imation scheme. In Section III, we describe the 2D realistic 
quantum dot model used for the study the scrambling and gate 
effects. The main results of this study are presented and ana- 
lyzed in Section IV. 



II. THEORY 

In both the KS-SDFT and S-TF approaches, for a system 
with N electrons and total spin S z , one solves a Schrodinger 
equation with a spin-dependent potential, 
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^(r)}C(r)=e^f(r) 



(1) 
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U cxt (r) + U H (v; [n]) + C/ x CT c (r; [n a , n 13 }) (2) 



where the total potential is the sum of the external, Hartree, 
and exchange-correlation contributions, and a = a, (5 de- 
notes spin-up and down, respectively. In KS-SDFT^ one 
knows that the potentials are functionals of the spin densities 
n a (r) — ^(r)! 2 , which are solved self-consistently 

under the constraint / n a (r)dr = N a with N a = (N + 2S Z ) /2 
and N 13 = (N — 2S z )/2. In analogy to the Koopmans's 
theorem in Hartree-Fock theory, it has been proved that the 
highest occupied Kohn-Sham orbital energy is identical to the 
chemical potential if the exact exchange-correlation potential 
is used. 20 This association provides a physical meaning to the 
self-consisent eigenvalues which can be contrasted with that 
of the reference Hamiltonian eigenvalues. 

In the S-TF approach, the basic idea is to start from 
a smooth semiclassical approximation, i.e. the Thomas- 
Fermi theory, and introduce quantum interference by consid- 
ering, first, single-particle corrections and, then, the effect of 
screened interactions between the oscillating part of the elec- 
tron density. In this case, then, one uses the TF potential in 
Eq. ([0. From the resulting ef and ip? one obtains an approxi- 
mation to the KS-SDFT total energy valid up through second- 
order in the oscillating density, n° sc (r) = n CT (r) — njp(r)iHii^ 
In this approximation, the relation between the KS-SDFT and 
Thomas Fermi single-particle levels can be expressed asi& 
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where the residual interaction terms 
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Se ni = , 



I drdr'|VA CT (r)|V s ^'(r,r')<: c (r') (4) 



corresponds to the interaction of ijjf with the density rip- 
ples due to interferences. The screened interaction potential 
V£.f (r, r') is expressed explicitly in terms of the functional 
derivatives of U a {v; [n Q ])>A& When for instance an electron is 
added into the quantum dot, both the genuine scrambling ef- 
fect (i.e. the variations of the Thomas Fermi levels) and the 
residual interaction terms Eq. (0} will affect the KS-spectrum. 

In this study we evaluate the scrambling and gate effects in 
both KS and S-TF approaches. For the S-TF case, the spin- 
dependence of the potential has little effect on its spectrum 
other than a constant shift, so we use the spin-independent 
Thomas-Fermi potential. For KS, we first solve the full spin- 
dependent KS equations, but calculate the scrambling and gate 
effects only from a orbital energies. Only minimal spin states, 
S z = for even N and S z = 1/2 for odd N, are considered. 
The spin indices will therefore be dropped in the remainder 
of the paper. The numerical methods that are used t o solve 
KS and TF equations are described in details in Refs l2ll22L 
respectively. 



III. MODEL SYSTEM 

The model system we use for investigating gate and scram- 
bling effects is a realistic 2D lateral quantum dot. 23 The elec- 
trons are at the heterointerface a distance Zd below the sur- 
face of the heterostructure. For the electrostatic potential, we 
use the mid-gap pinning model for the boundary condition at 
GaAs free surface^i 2 ^ 2 ^ We impose Dirichlet boundary con- 
ditions on the top surface and Neumann conditions in infinity 
in all other directions, 23,27 allowing the external potential to 
be calculated from 
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(5) 

where f7 g (r) is the electrostatic potential on the top gate 
surface, 26 and Uqw(z) is the confining potential in the growth 
direction due to the quantum well structure from which the 
quantum dot is fabricated. In addition, the Hartree poten- 
tial has an image term due to the coupling with the top 
surfaceriST. 



U H (r; [»]) = jdr'n(v') 



1 



(|r - r'| 2 + 4^2)1/2- 



(6) 

For a complete description of our treatment of the electrostatic 
potential see Ref. [23[ 

The shape of the top confining gate used here is shown in 
Fig|2j it is designed to model typical irregular dots investi- 
gated experimentally. 6 Negative voltages are imposed on the 
top (Z7t), bottom (Ub), and plunger (U p ) electrode gates; the 
electron number in the dot is controlled by U p . The single- 
particle dynamics of the system is expected to be chaotic, 
which is confirmed by the agreement between the nearest- 
neighboring spacing distribution of the single-particle levels 
and the Wigner surmise distribution! 2 ^ 

We now introduce quantities to characterize the scrambling 
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FIG. 2: (Color online) Left panel: Schematic of the shape of the top 
confining gate used in this study. Negative voltages are imposed on 
the shaded region, depleting the electrons underneath so that the mo- 
tion of electrons is confined to a small region. To obtain enough 
statistics, 24 different irregular confining shapes are generated by 
taking different values of Xi, X2, Y\ and Y2. Right panel: The near- 
est neighbor spacing distribution of the single-particle levels calcu- 
lated in the Thomas-Fermi effective potential (histogram) compared 
to the Wigner surmise distribution (line). 

and gate effects. For the scrambling effect, we define 

5ei(SN) = £l (N^ + SN) - e. ( (iV (0) ). (7) 

where £i(N) is the i lh single-particle energy in either the TF or 
KS effective potential with N electrons. Scrambling is quan- 
tified by the magnitude of fluctuations in Sei(SN): 

* sm ^«{ S ^-^ SN »] (8) 

where (Sei(5N)) denotes a linear fit of S£i(8N) as a func- 
tion of i, and A is the mean level spacing. The gate effect is 
similarly characterized as 

a g (SU p )^a{ S£ ^-^ SU ^) (9) 

with Se^SUp) = £l (C/j 0) + SU P ) - e l (U^ )) ), and {Se^SUp)) 
its linear fitting. It is more convenient to write a g as a func- 
tion of SN* = SUp/ (SU°), where (5U°\ is the average con- 
ductance peak spacing. SN* can be regarded as the induced 
electron number due to a change of the gate voltage. 

IV. RESULTS AND DISCUSSION 

The scrambling and gate effects are mixed together in CB 
conductance peak spacings. It is desirable, however, to first 
study them separately; while the separation of the two effects 
is difficult to implement experimentally and requires sophis- 
ticated design, 6 it is straightforward in numerical investiga- 
tions. For the scrambling effect, we fix the external confining 
potential and calculate the TF and KS single-particle spectra 
at each N e [50, 70]. For the gate effect, we fix N = 70 and 
S z = 0, and scan the plunger gate voltage for SN* up to about 
20. Statistics are obtained by averaging over different levels 
and 24 confining gate shapes. 

Fig- 0a) shows the scrambling effect calculated from both 
TF and KS single-particle spectra. A remarkable feature is 
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FIG. 3: (Color online) Scrambling magnitude as a function of SN 
(upper) and gate fluctuation as a function of 57V* (lower) from the 
TF (filled) or KS (empty) single particle spectra with z& = 15 nm 
(circles) or 50 nm (triangles). Inset: The result of a universal random 
matrix model. 

the large difference between TF and KS results. a s from TF 
spectra increases smoothly as a function of SN, but is smaller 
than 0.15 even at SN = 20. 29 In contrast, a s from KS spectra 
increases rapidly at first, and saturates for SN > 4. For SN > 
8, <7 S shows some modulations, the reason for which is not yet 
clear. The effect of the dot depth, z& = 15 nm vs. 50 nm, is 
quite small. The fact that the KS results are much larger than 
TF ones shows that the residual interaction effects included in 
the self-consistent KS energies [Eq. ] dominate the genuine 
scrambling evaluated with TF. 

Fig-EJb) shows the gate effect calculated from both TF and 
KS single-particle spectra. The gate effect in TF and KS spec- 
tra are qualitatively the same; in both cases a g first increases 
and then saturates at some large 5N* (about 10 for z c \ = 15 
nm in the TF case). The saturated values of a g are about 0.6, 
larger than that of c s (KS). The gate effect is quite sensitive to 
the depth of the dot, especially for the TF case; it is larger for 
shallow dots, as expected because the gate becomes sharper 
and better defined. Finally, note that the gate effect is about 
one order of magnitude larger than the TF scrambling effect. 

To further understand these results, we model the paramet- 
ric evolution of the single-particle spectrum by using a ran- 
dom matrix Gaussian processpS 

H{x) = cos(xir/2)H 1 + sm(xTT/2)H 2 (10) 

where Hi and H2 are independent random matrices belong- 
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ing to the Gaussian orthogonal ensemble (GOE). 2,28 We define 
a u (x) similarly to Eqs. or to characterize the change 
in the single -particle spectrum of H (x) due to the variation of 
x. The inset of Fig. 13b) shows a u vs. x obtained from 500 
realizations of 30 x 30 random matrices. a u saturates at about 
x = 0.4, and the saturated value is about 0.6 in agreement with 
that of er g . We notice that the saturation behavior of a s from 
KS levels is very different from that of er u ; this is because the 
KS single-particle levels contain implicitly some interaction 
effects so that their variation with N is not a simple Gaussian 
process. 

The agreement between the KS and TF results for the gate 
effect in Fig. |3jb) is particularly striking when compared to 
the sharply different magnitudes for the scrambling effect in 
Fig- 0a). The underlying reason for this difference between 
the gate and scrambling effects is the extra electron added in 
the case of scrambling: First, the contribution to the fluctua- 
tions of by most of the residual interactions is small be- 
cause of the small change in wave functions caused by the 
additional potential. Second, the contribution by residual in- 
teractions involving the added electron j are, however, much 
larger: the fluctuations here involve the deviation of |-0J(r)| 2 
from the smooth density rather than the small change in an 
already filled level. Thus the fluctuation of the residual inter- 
action contribution to e^f is substantially stronger in the case 
of scrambling than for the gate effect. 

We now turn to actual spacings between CB conductance 
peaks and calculate the scrambling and gate effects. The po- 
sition for the N th peak, Ug, at which the energies for N — 1 
and N electrons are equal, is determined by n(N,U g ) = 
E KS {N, U g ) - E KS (N - 1, U g ) = if chemical potentials 
in the leads are taken to be zero. From the TF spectrum at 
this U g , the scrambling and gate effects are the standard devi- 
ation of e t (N, uf ] ) — Si(N—l, U { g N) ) and e t (N, U [ g N+1) ) - 
£i (N, Ug ), respectively. For Zd — 15 nm, the former is equal 
to 0.009A, and the latter 0.07A, using the same confining gate 
shapes and parameter ranges as above. Note the good agree- 
ment with the results in Fig.|3]for SN = 1 or SN* = 1. 

Comparison with earlier evaluations of these effects yields 
important insights. We start with the scrambling effect ; in 
particular, using the expressions derived in Refs. Il3ll4l for 
scrambling associated with the Thomas-Fermi spectrum - 
which we think of as the "genuine"scrambling effect - yields 
af red (<W = 1) ~ 0.06 for a dot with N = 70 electrons. This 
prediction is six times larger than the value obtained here. 
There are two main differences between the earlier situation 
and ours: the confining potential here is smooth while it was 
assumed to be hard-wall in Refs. Il3ll4l and we effectively 
have a top gate across the whole dot because of our boundary 
condition. The insensitivity of the scrambling magnitude to 
the spacing Zd suggests that the top gate has little effect. We 
therefore conclude that there is significantly less scrambling 
in a smooth confining potential than in a hard wall dot. 

Further support for this conclusion comes from the change 
in potential upon adding an electron: the hard-wall gives rise 
to a square-root singularity in this quantity; 11 the absence of 
a singularity in the smooth case naturally leads to a weaker 



effect. As a consequence, the scrambling in the KS spectra, 
which contains both genuine scrambling and some residual 
interaction fluctuations [see Eq. (|3}], appears to be entirely 
dominated by the latter. Evaluating the magnitude of these 
residual interaction fluctuations as an extra electron is added, 
following the semiclassical random plane wave approach in 
Refs. ll3ll4l gives 0. 18 A for N = 70. Taking into account that 
one is not very far into the semiclassical regime, this is quite 
compatible with the value a s (8N* = 1) » 0.13 computed for 
the Kohn-Sham spectra [Fig.[3ja)]. 

For diffusive transport in a weak disordered potential, 
scrambling has been studied using the statistics of single- 
particle wavefunctions in that case. 2 In that context, scram- 
bling grows linearly with SN while residual interaction ef- 
fects grow as (SN) 1 / 2 . It is interesting to note that our data in 
Fig. |3j a) for a ballistic dot show roughly the same behavior. 

With regard to the gate effect, it was argued that a "generic" 
gate should have the same effect as TF scrambling 14 . For the 
gate considered here, we however observe a significant differ- 
ence. A "generic" gate is one which couples approximately 
uniformly to the dot - a back gate, for instance. The oppo- 
site extreme is a gate coupled very locally to a point in the 
dot, thus producing a rank one perturbation. Such perturba- 
tions are known to completely decorrelate the spectra for the 
phase shift it/2 necessary to add an extra electron. 31 Here we 
see that lateral plunger gates are an intermediate case: they 
produce fluctuations which are significantly stronger than the 
scrambling effect but remain moderate on the scale of the 
mean level spacing. 

The gate effect here is nevertheless strong enough to pro- 
duce spin transitions. Recently Kogan, et al. reported ex- 
perimentally a singlet-triplet transition in zero magnetic field 
driven by changing the confining potential. 32 In our calcula- 
tions, we also observe spin transitions caused by gate-induced 
shape deformation. The probability of a transition - the frac- 
tion of peak spacings in which the ground state spin flips as 
Up changes - is 5.8% for Zd = 15 nm and 4.6% for z<j = 
50 nm. In addition to the simple singlet-triplet transition, dif- 
ferent patterns appeared in our calculations. In many cases, 
the spin flips are paired: in a single spacing interval, the spin 
changes from one value to another and then back again, which 
agrees with the general picture connecting spin transitions 
to avoided-crossings caused by shape-deformation. 33 In some 
rare cases, triple-transitions are observed, such as S z = 2 — ► 1 
->0->l or 3/2 ->l/2 ->3/2 -vl/2. In other cases, the spin 
transition is unpaired, presumably because pairing due to an 
avoided-crossing is destroyed by a change in electron number. 

Several experiments have traced the correlation between 
ground states or excited states as the number of electrons 
changespiiMi and often see a surprising degree of correla- 
tion. The relatively weak scrambling and gate effects that we 
find offer a way to understand these results. 

To conclude, we have investigated the scrambling and gate 
effects for two different one-particle spectra - self-consistent 
Kohn-Sham and Thomas-Fermi - of a realistic model quan- 
tum dot. Our main findings are: (1) The genuine scrambling 
effect - the one associated with the Thomas-Fermi spectra - is 
significantly smaller for the smooth potential considered here 
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than in earlier work using hard wall confinement. (2) As a 
consequence, scrambling for the Kohn-Sham spectra, which 
involves both genuine scrambling and residual interactions, 
is entirely dominated by residual interactions. Its magnitude 
agrees with those from a random plane wave model of the 
wave functions. (3) Finally, fluctuations caused by the gate 
are similar in magnitude for the two spectra. A lateral plunger 
gate causes significantly larger fluctuations than those caused 
by scrambling, contrary to the case of a uniform "generic" 



gate. The magnitude of the gate effect remains nevertheless 
moderate compared to the mean level spacing. 
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